arXiv: 1504.04268v2 [hep-ph] 19 Jan 2016 


Magnetized plasminos in cold and hot QED plasmas 


N. Sadooghl!] and F. Taghinava^ 

Department of Physics, Sharif University of Technology, P.O. Box 11155-9161, Tehran-Iran 

The complete quasi-particle spectrum of a magnetized electromagnetic plasma is systematically 
explored at zero and nonzero temperatures. To this purpose, the general structure of the one-loop 
corrected propagator of magnetized fermions is determined, and the dispersion relations arising from 
the pole of this propagator are numerically solved. It turns out that in the lowest Landau level, 
where only one spin direction is allowed, the spectrum consists of one positively (negatively) charged 
fermionic mode with positive (negative) spin. In contrast, in higher Landau levels, as an indirect 
consequence of the double spin degeneracy of fermions, the spectrum consists of two massless col¬ 
lective modes with left- and right-chiralities. The mechanism through which these new collective 
excitations are created in a uniform magnetic field is similar to the production mechanism of dy¬ 
namical holes (plasminos) at finite temperature and zero magnetic fields. Whereas cold magnetized 
plasminos appear for moderate magnetic fields and for all positive momenta of propagating fermions, 
hot magnetized plasminos appear only in the limit of weak magnetic fields and soft momenta. 

PACS numbers: ll.lO.Wx, 11.30.Na, 12.38.-t, 12.38.Aw, 13.40.-f 


I. INTRODUCTION 

Research on matter under extreme conditions has 
provided new insight in the physics of heavy-ion colli¬ 
sions (HIC) and the astrophysics of compact stars. Ex¬ 
treme conditions consist of high temperature, large den¬ 
sity and/or the presence of intense magnetic fields. The 
latter affects, in particular, the phase diagram of Quan¬ 
tum Chromodynamics (QCD), and plays a significant 
role in the dynamics of relativistic fermions at zero and 
nonzero temperatures (for recent reviews, see Qi). The 
phenomena driven by external magnetic fields have also 
various applications in quasi-relativistic condensed mat¬ 
ter systems, such as graphene and Dirac semi-metals (see 
3 and the literature therein). 

Most theoretical studies deal with the idealized limit 
of constant and homogeneous magnetic fields. Standard 
field theoretical methods thus lead to the exact solution of 
the relativistic Dirac equation in a uniform magnetic field. 
One of the main consequences of the presence of constant 
magnetic fields is a certain dimensional reduction of the 
dynamics of propagating fermions in the lowest Landau 
level (LLL). The latter leads to a dynamical mass gener¬ 
ation, and enhances the production of chiral condensates. 
This phenomenon, which is known as magnetic cataly¬ 
sis ii, modifies, in particular, the phase diagram of 
QCD in the chiral and color superconductivity phase Q. 
Another effect is the appearance of certain anisotropies 
in the dynamics of magnetized fermions in the longitu¬ 
dinal and transverse directions with respect to the direc¬ 
tion of the external magnetic field. These anisotropies 
include those in the neutrino emission from magnetars 
0, or anisotropies arising in the group velocities, refrac¬ 
tion indices and decay constants of mesons in hot and 
magnetized quark matter Q. Recently, the anisotropy 
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appearing in the equation of states of magnetized quark 
matter is studied in 

A similar privileged reference frame is also defined by 
a heat bath, and affects, in particular, the quasi-particle 
spectrum of electromagnetic and quark-gluon plasmas at 
finite temperature. Non-trivial bosonic and fermionic col¬ 
lective excitations, such as plasmons and plasminos are 
shown to be dynamically generated in hot QED and QCD 
plasmas, in addition to the normal bosonic and fermionic 
modes [lOlllll- In particular, plasminos are known to be 
collective excitations that arise as one of the poles of the 
one-loop corrected fermion propagator at finite tempera¬ 
ture [12| . In the chiral limit, plasminos are characterized 
by their negative helicity to chirality ratio, in opposite 
to that of normal modes. They are intensively studied 
in the context of Yukawa theory, QED and QCD [T3l - [l^ . 
In 1^, for instance, it is shown that the contributions 
from plasminos modify the transport properties of rel¬ 
ativistic plasmas in- and out-of equilibrium. They also 
lead to the appearance of sharp structures (singularities 
and gaps) in the decay [l^ and production rates [l^ of 
particles produced in relativistic and ultra-relativistic col¬ 
lisions. These structures provide unique signatures for the 
presence of deconfined collective quarks in the plasma of 
quarks and gluons. 

In principle, nontrivial collective modes can also be 
created in magnetized plasmas through the same mech¬ 
anism as the one leading to the appearance of collective 
modes at finite temperature and zero magnetic fields. It 
is the purpose of the present paper to look for possi¬ 
ble dynamical generation of fermionic excitations (plas¬ 
minos) in electromagnetic plasmas in the presence of con¬ 
stant magnetic fields, and study their properties at zero 
and nonzero temperatures. Apart from various other ap¬ 
plications, plasminos may play an important role in the 
physics of HICs. Very strong magnetic fields, which, ac¬ 
cording to recent experimental results, are believed to be 
created in early stage of non-central HICs da may af¬ 
fect, among others, the energy dispersion of deconfined 
quarks. The latter are believed to be produced in the 




quark-gluon plasma in the same stage as the magnetic 
fields. It is therefore important to explore the quasi¬ 
particle spectrum of Dirac equation in the presence of 
external magnetic fields at finite temperature, and study 
the properties of the potentially created collective modes 
under these conditions. 

In this paper, we will particularly focus on the mech¬ 
anism of the production of plasminos in the presence of 
external magnetic fields. We will, in particular, deter¬ 
mine the general structure of the one-loop corrected prop¬ 
agator of magnetized fermions, and, following the histor¬ 
ical path that has led to plasminos at finite temperature 
[l^ , solve the dispersion relations arising from the pole of 
this propagator. To this purpose, we will first determine 
the general structure of the tree-level fermion propaga¬ 
tor in the presence of a constant magnetic field by mak¬ 
ing use of the Ritus eigenfunction method [3. The free 
fermion propagator will then be combined with the one- 
loop self-energy of magnetized fermions. This will result 
in the desired one-loop corrected propagator of magne¬ 
tized fermions. We will show that, in contrast to LLL, 
the dressed fermion propagator in higher Landau levels 
(HLL) can be decomposed into two parts, each of them 
leading to a separate energy dispersion relation for mag¬ 
netized fermions. This fact, which eventually leads to the 
appearance of magnetized plasminos, is an indirect con¬ 
sequence of the double spin degeneracy in HLL, in con¬ 
trast to LLL, which is occupied with only one positive or 
negative fermion with positive or negative spin. We will 
show that at finite temperature and in the limit of weak 
magnetic fields, where HLLs have also to be taken into 
account, the spectrum consists of two massless collective 
modes with left- and right-chiralities. Moreover, it can be 
shown that whereas cold magnetized plasminos appear for 
all positive momenta of propagating fermions and moder¬ 
ate magnetic fields, hot magnetized plasminos appear only 
in the limit of soft momenta and weak magnetic fields. 

The organization of this paper is a follows: In Sect, 
im we will review two independent topics related to the 
main subject of the paper: In Sec. Ill A[ we will first show 
how thermal plasminos arise from the pole of the one-loop 
corrected fermion propagator in the Hard-Thermal Loop 
(HTL) approximation [l9|. In Sec. IIIBI we will then 
review the Ritus eigenfunction method [1^, and present 
the general structure of the free propagator of magnetized 
fermions in the momentum space. The general structure 
of one-loop fermion self-energy at zero and nonzero tem¬ 
peratures will be derived in Secs. IIII Al and IHIBI respec¬ 
tively. In Sec. IIII HI in particular, a certain HTL approx¬ 
imation in a constant magnetic field will be introduced, 
and the one-loop fermion self-energy will be presented in 
terms of a number of coefficients up to some integrations 
and a summation over Landau levels. In Sec. ITVl the 
general structure of the one-loop corrected propagator of 
magnetized fermions will be determined by combining the 
free fermion propagator from Sec. IIIBI and the one-loop 
fermion self-energy from Sec. IHII Here, two different cases 
of massive fermions in LLL and massless fermions in HLL 
will be considered, and the properties of the fermionic 


excitations in these two cases will be systematically stud¬ 
ied. In Sec. El after numerically determining the afore¬ 
mentioned coefficients, appearing in the one-loop fermion 
self-energy, we will study the spectrum of the fermionic 
excitations at finite [Sec. IV A| and zero [Sec. IV B| tem¬ 
peratures. Section ED is devoted to a brief summary of 
our results and a number of concluding remarks. 

II. REVIEW MATERIAL 

The main goal of the present paper is to determine the 
spectrum of Dirac particles in a constant magnetic field 
at zero and nonzero temperatures. To this purpose, and 
in order to fix our notations, we will briefly review, in 
this section, two independent topics related to the main 
subject of this paper. In Se c. Ill Al we will first repeat the 
computation presented in 07 and introduce the plas- 
mino excitations in a hot QED plasma. To do this, we 
will compute the one-loop correction to the fermion self¬ 
energy at finite temperature in a HTL approximation , 
and eventually determine the Dirac spectrum at finite 
temperature by solving the corresponding energy disper¬ 
sion relation. In Sec. Ill HI we will then briefly review the 
Ritus eigenfunction method [l^ , and after presenting the 
free propagator of magnetized fermions in the coordinate 
space, will determine its general structure in the momen¬ 
tum space. Our review will mainly base on the Ritus 
eigenfunction method, presented in (2^, and generalized 
in Q for a multi-flavor system of charged fermions. 

A. Plasminos in a hot QED plasma 

In [T^ . the Dirac spectrum is determined for a hot 
QCD plasma. In this paper, however, we will focus, for 
simplicity, on the C/(l) subgroup of the SU{N) gauge 
group. The difference between our 17(1) and the SU{N) 
case, discussed in [HIS, is a certain prefactor Cp, the 
quadratic Casimir constant of the fermion representation, 
in the corresponding expression to the one-loop fermion 
self-energy. For the Abelian U{1) gauge group, (7^ = 1 
and for the fundamental representation of SU{N) gauge 

N ^-1 

group, Cf = 


k-p 
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FIG. 1: One-loop fermion self-energy 

Let us start by determining the one-loop fermion self¬ 
energy of massless fermions at finite temperature.^ Hav- 


^ At high enough temperature, the mass of fermions can be ne- 
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ing in mind that the interaction term of fermions and pho¬ 
tons in QED is given by = —e'07^^^'0, the fermion 
self-energy E(fc) at finite temperature is given by [see Fig. 

H] 


+ 00 


E(fc) = e^T 

n— — oo '' 
+00 
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(II.l) 


Here, Df^^(k) = ^ and S{p) = ^ are free photon 
and fermion propagators, respectively.^ In the imagi¬ 
nary time formalism, po is to be replaced by iujm where 
ojn = (2n + 1)ttT is the fermionic Matsubara frequency. 
Since in the presence of a hot medium, the relativistic 
Lorentz invariance is broken by introducing the reference 
frame corresponding to the heat bath, it is necessary to 
perform the following separation between the temporal 
and spacial components of T,{k) 


E{k) = 7 oSt(fc) - 7 • Ss(fc). 
Introducing Ei = jpj, F 2 = |p — k|, and using 

Af{po,Ei) = 


(II.2) 


Ab{ko-po,E2) = 


Pl + Er 

1 


(II.3) 


{ko - PoY + El ’ 

the temporal and spacial components of E(fc) are given 

by 

(Pp 


Etik) = - 2 e^y 

(Pp 


Yls{k) = - 2 e^ 


(27r)3 


P Is{k,v), 


(11.4) 


with 


It = T ^ iuJnAf{iuJn,Ei)Ab{ko - iuJn,E2), 

n— — oo 
00 

Is = T A f{iuJn,Ei)Ab{ko - iuJn,E 2 ). (II.5) 


n— — oo 


Following the standard method introduced in [2l|, the 
summation over Matsubara frequencies can be performed 
by making use of 

52 l + fbis2E2)-ff{siEi) 

4:E2 ko — siEi — S2E2 

S 1 S 2 1 + fb{s2E2) - ff[siEi) 
4:EiE 2 ko — siEi — S2E2 


it(fc,p) = - 


Si ,S 2 =± 


Is{k,p) = - 


Si,S2=± 


(II. 6 ) 


Here, the fermionic and bosonic distribution functions are 
defined by 

= e/3^0 + 1 ’ 

fb{io) = ^ 1 (m) 


with P = T“^, and £0 = i|.^| for massless fermions. Re¬ 
placing the above expressions in (III. 41) . performing an ap¬ 
propriate expansion in jkj ^ jpj with jpj ^ Tp and keep¬ 
ing only the leading contributions from the resulting 
expressions, St and Sg in this HTL approximation are 
given by 


St(fc) 

Ss(fc) 


2 |ki Vfco-iki;’ 

1 Ao + 

|k| V 2 |k| \ko-\k\)) 


k. 


(H. 8 ) 
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Here, mj^ = is the Debye mass. To perform the inte¬ 
gration over p (here denoted by p), the following integrals 
are used 



dppffip) = 


dppfbip) = 


6 ’ 
12 ■ 


(H.9) 


Plugging at this stage the results from (jH.811 in (|H.2p , and 
combining the resulting expression for S(/c) with the in¬ 
verse of the free fermion propagator, S~^{k) = 70 ^ 9 — 7 -k, 
the inverse of the one-loop corrected fermion propagator, 
5-1 (fc) = S-^{k) - E{k) reads 


S = 7 oHo - 7 • k Hg, (11.10) 


with 


m 




|k| 


ko - |k|/ ’ 




ko - Ik| 


( 11 . 11 ) 


The dressed fermion propagator up to one-loop order is 
therefore given by 


S{k) = 


70 - 7 • k 7o -I- 7 • k 


2£»4 


2 £»_ 


( 11 . 12 ) 


with 


D± = AotAs, (11.13) 

and Ao/s from (III.111) . To determine the spectrum of 
Dirac particles, we set either D± = 0 or det[iS“^(fc)] = 0, 


3 In this approximation E 2 is, in particular, replaced by i?2 ~ |p| ~ 
|k| cos 0 , with 0 the angle between k and p. 


glected. 

^ In what follows, we will denote the free and dressed fermion prop¬ 
agators with S and <5, respectively. 
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and arrive at two solutions Aq = or more explicit¬ 
ly at two different energy branches for particles (p) and 
dynamical holes (h) [T^ . 


Er. = |k| 


m 

Tk 


^ 1 + I (1 - In 


Zp + 1 

Zr,- 1 


Eh = 






(11.14) 


where Zp/h = Replacing the expression on the 

right hand side (r.h.s.) of Eh with 


®‘ = |k|““'(k' + ikr^)’ 

the dependence of dimensionless quantities ko/mo for 
ko = Ep and fco = Eh on |k|/TOD can be determined 
numerically. Two different energy branches correspond¬ 
ing to D+ and Zl_ arise. They are demonstrated in Fig. 

12 
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FIG. 2: (color online). Two different energy branches, corre¬ 
sponding to the energy dispersion relations D+ = 0 of par¬ 
ticles (blue curve) and D- = 0 of holes (red curve). The 
denominators D± are defined in (III.1311 . The particles/holes, 
whose energy dispersion relations are demonstrated by these 
two branches, have positive/negative helicity {'H) to chirality 
(x) ratio. 


In what follows, we will argue that the energy branch 
arising from D+ = 0 {D- = 0) [blue (red) curve in Fig. 
[2] corresponds to a particle (hole) with positive (nega¬ 
tive) helicity (H) to chirality (x) ratio. To do this, let us 
determine the eigenvectors of the numerators in (III.12P , 


A4 = 7o T 7 • k = 


0 1 T rr • k 

l±(T-k 0 


(11.15) 


Here, cr = (CTi,CT 2 ,cr 3 ) are the three Pauli matrices.® For 
branch, A/+ has the following two nontrivial eigenvec¬ 
tors 


= 

= 

Defining the helicity 


ki + ik2 
ks - |k| 
ki + ik2 


, 1 , 0,0 


(11.16) 


H = S-k, (11.17) 


with S = diag ((T,cr), and right- (i?) as well as left- (L) 
chirality operators 



1 L- 2 ’ 

(H.I8) 

ii) 

it turns out that W_)_ , f = 

1 , 2 satisfy 



= +wi^\ 



= - wf. 

(H.I9) 


Having in mind that right (left)-handed particles have 
positive (negative) chirality, x = +1 (x = ~1)) the re¬ 
lations in (III.191) show that fermions, whose energy dis¬ 
persion relation is given by branch in Fig. [2l have 
a positive helicity to chirality ratio, E = _|_i_ They are 
therefore particles. Similarly, it can be shown that for 
D-, the eigenvectors of Af- are given by 


= ( 

, ’ ^kl+ik2’ ) ’ 

= ( 

+ tool 

7 l+*fc2’ ’ ’ )' 

which satisfy 


VRW^y = w^y, 

nwy = -wy, 

Vhwy = wy, 

nwy=+wy 


( 11 . 20 ) 


( 11 . 21 ) 


These relations indicate that fermionic mode, whose en¬ 
ergy dispersion relation is given by Zl_ branch in Fig. [5J 
have a negative helicity to chirality ratio, ^ = —1. Ac¬ 
cording to the standard terminology of relativistic quan¬ 
tum mechanics, they are therefore holes. In conclusion, 
the spectrum of Dirac particles at high enough tempera¬ 
ture, consists of two different soft excitations with positive 
and negative helicity to chirality ratio. In Sec. IHII the 
above HTL approximation will be used to determine the 
general structure of one-loop self-energy of fermions in 
the limit of soft momenta and weak magnetic fields. This 
will eventually lead to hot magnetized plasminos in these 
limits [see Sec. E]- 


^ To derive iUltsIi . we have chosen the following chiral representa¬ 
tion of Dirac 7 -matrices 


7^ = 


/ 0 (je 'x 

5 / 

-10 7 

V se 0 ; ’ 

7 ®=( 

, 0 1 ) 


^ Here, fco is replaced by Ep [first expression in and Eh 

[second expression in fiTTili ]. 
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with o-e = o") and = ( 1 , —cr). 
































B. Magnetized fermions in a cold QED plasma; The 
Ritus eigenfunction method 

In this section, we will solve the Dirac equation of 
positively and negatively charged fermions in the presence 
of an external magnetic field 

(7 • - mg ) i ) = 0, (11.22) 

with = idfj, + eqA‘^^, and rriq the fermionic mass. 
Here, e > 0 is the unit electric charge, and q = ±1 indi¬ 
cate the positive and negative charges of Dirac fermions. 
The gauge field = (0, 0,Ha:i,0) is chosen so that it 
yields a magnetic field B, aligned in the positive third di¬ 
rection, B = Be^ with B > 0. In the method originally 
introduced by V. I. Ritus in [T^ , and generalized recently 
in [1, [23 for multi-flavor systems, (III.22|) can be solved 
by making use of the Ansatz for a Dirac 

fermion with charge q. The Ritus eigenfunction then 

satisfies 

(7.n)E(«) =EW(7.pp), (11.23) 

and the free Dirac spinor up satisfies (^p — m)up = 0. As 
it turns out, the Ritus momentum pp is given by 

Pp = {po,0,-Sq^/2p\qeB\,p3), (11.24) 


iq = {\qeB\) and the normalization factor Up = 
{2'Pp\^lq) In (III.301) . Hp{x) is the Hermite poly¬ 
nomial of order p. Using the above results, it is easy to 
determine the free propagator of Dirac fermions in the 
coordinate space at zero temperature [1, , 



(II. 3 I) 


Here, Pp and Pp'^'^ for Dirac fermions with charge q are 
given in (III.241) and (IH.26|) . respectively. 

In what follows, we will derive the general structure of 
free fermion propagator in the presence of a constant mag¬ 
netic field at zero temperature in the momentum space. 
To do this, let us first consider the free fermion propa¬ 
gator (III.311) in the coordinate space. The free fermion 
propagator in the momentum space is given by perform¬ 
ing an appropriate inverse Fourier transformation, where, 
instead of the standard plane wave basis, the Ritus basis 
(|H.25|) is used 


(fc, k') = I d^xd^y E(«) {x, (x, y)E(f)t (y, k'). 

(11.32) 


where p labels the Landau levels in the external magnetic 
field B, and Sq = sgn{qeB). The Ritus eigenfunction for 
Dirac fermions with charge q can be derived from (III.231) 
and (jH.24() . and reads 

E^;\x,p) = e-*P'®P^«)(xi,p2), (11.25) 

where p= {po,0,P2,P3), x= (xo, 0, X 2 , X 3 ), and 

pig) = pD) f+s, ^ , (11.26) 

with = V± and as well as 


V± = 


I ± 17172 
2 


(11.27) 


In (III.261) . Hp = 1 — i5po considers the spin degeneracy 
in the Landau level p (see below). For eB > 0, Pp^^ 
for positively {q = -fl) and negatively (9 = — I) charged 
particles read 


= ^±/p^ + npP^/pT, ( 11 . 28 ) 


respectively. Here, f^{xi,p 2 ) are given by 

/+ = $p(xi-£2p2), for p = 0,1,2,---, 

f- = $p_i(xi-^2p2), for p = 1,2,---, 

(11.29) 

with 

4>p(x) = Op exp Hp (^y) ’ 


Plugging En\x,k) from (IH.26I) on the r.h.s. of (IH.32|) . 
and performing the integrations over x = (xq, 0, X 2 , X 3 ) as 
well as y = (po, 0 , j/ 2 , 2/3), we arrive first at 

S'W(/c,fc') = {2T:f5^{k - k')S‘'^\k), (11.33) 

with 

00 . 

Ck) = E ^Pn ■ (11-34) 

p=0 rp ^9 

Here, Wpn is defined by 



Using at this stage the definition of Pp‘^^ from (IH.26p . 
the properties of the projectors = P^'^ and 

pD)pD) _ integration 



dxifp‘'‘‘[xi,k2)ft:‘'’‘{xi,k2) = Spn, 


(11.36) 


which arises from the standard integration over Hermite 
polynomials 


P + CXD 


dx e ^ Hp{x)Hn{x) = 2Pp!v^5p„, (11.37) 


we arrive, after some calculation, at 

= [pI"^ + npH^pi^^pn. (H.38) 
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Plugging then (IIL38I) in (III.31) . splitting into longitudinal 
and transverse parts, ^ with = 7II • ky 

and ) 4 _|_ = 7'*‘ • kj^, and ky = (fco,0,0, fca) and kj_ = 
(0,0,-SqY^2p|ge.B|,0) as well as 7 y = ( 70 , 0 , 0 , 73 ) and 
7 j_ = ( 0 , 71 , 72 , 0 ), and using eventually® 

pOj) 7^11 = 7/^11 ( 1139 ) 

the free fermion propagator (|II.3ip in the momentum 
space reads 




(n + + mg )+ 

kl — mi 


(11.40) 


Here, H^ = n„ and = 1 are also used. Accord¬ 

ing to (|II.40I) . the free fermion propagator for positively 
and negatively charged particles in the LLL and HLL at 
zero temperature reads 


For n = 0 
For n ^ 0 


Si^\k) 


_ 'P± 

1 ^11-m, ’ 
1 

kn — 'mq ’ 


(11.41) 


where the superscripts ± on stand for q = ±1. The 
appearance of the projectors V± in Sq reflects the fact 
that in the LLL, the charged fermions have either a posi¬ 
tive (positively charged fermions) or a negative spin (neg¬ 
atively charged fermions), while in HLL fermions a certain 
double spin degeneracy occurs. In other words, in HLL 
positively (negatively) charged particles have both pos¬ 
itive and negative spin orientations. Let us also notice 
that, the appearance of )4y in the denominator of 
is related with the expected dimensional reduction from 
Z1 = 4 to H = 2 in the LLL Q. In the next section, we 
will derive the general structure of the one-loop perturba¬ 
tive correction to the fermion self-energy of positively and 
negatively charged particles for H ^ 0 and at T ^ 0. The 
results will then be combined with (IH.41I) , and lead to the 
general structure of the one-loop corrected propagator of 
magnetized fermions at zero and nonzero temperatures. 


III. ONE-LOOP SELF-ENERGY OF 
MAGNETIZED FERMIONS AT ZERO AND 
NONZERO TEMPERATURES 


In this section, we will determine the one-loop self¬ 
energy contribution of charged fermions in the presence 
of a constant magnetic field at zero and nonzero temper¬ 
atures. In Sec. IIII Al we will first consider the zero tem¬ 
perature case, and will show that the one-loop fermion 
self-energy, 'En\k), has the following general structure 

s(9) = ^y(P+A^«^ 

+mg{r+d_^^ (iii.i) 


In llll.39t . /r|| and denote /i = 0, 3 and fj, = 1,2, respectively. 


where the superscripts (g) stand for positively {q = -|-I) 
and negatively {q = —I) charged fermions. According to 
the notations introduced in the previous section, ky = 

(fco,0,0, fcs) and kj^ = (0, 0, -Sq-\/2n|geB|, 0). 

In Sec. IhTrI we will then show that at finite tem¬ 
perature an additional splitting between the zero and the 
third components of ky occurs. The general structure of 

'sl?\k) in a hot and magnetized QED plasma therefore 
reads 

+M'P+B^+ 

(HL2) 

Let us notice that the additional splitting between the 
zero and third components occurs because of broken 
Lorentz invariance induced by the heat bath. 

In this section, the coefficients and from 

Hum and and from (THL^ will be 

analytically determined up to a summation over Landau 
levels and a number of integrations. The latter will be 
then numerically performed in Secs. IV Al and IVHI where 
the spectrum of Dirac fermions in a magnetized QED 
plasma will be determined at zero and nonzero tempera¬ 
tures. The results will eventually be compared with the 
results presented in Sec. Ill Al for hot QED plasma and for 
vanishing magnetic fields. 


A. One-loop self-energy of magnetized fermions at 
zero temperature 


The one-loop self-energy of magnetized fermions at 
zero temperature is given by combining the free fermion 
propagator (IH.3I|) . the free photon propagator = 

in the Feynman gauge and the corresponding vertex 
—eq'y^ of a photon and a magnetized fermion pair. 


1=0 


dp dp 

(27r)3 (27r)4 


X {xi , g 2 ) g ^ 'Y' Pi^' {yi ,P2)7" 

\Pl '"'q) P 


+ rnq) 




/ 


(HL3) 


To derive the fermion self-energy in the momentum space, 
we use, as in the previous section, an appropriate inverse 
Fourier transformation, using the Ritus eigenfunctions ]E„ 
from (|H.25I) . 

= J d'^xd'^y {x, {x, (y, k'). (HL4) 


After performing the integration over x and y, we arrive 
at 


EW(fc,fc') = {2Tif 5\k - k'^Hk), (III.5) 
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with 


and 


= eS‘Y. / 7 SI 7 

f =0 ^ 

_ {Xi,yi]k2,P2) _ 

^ [(k|| - P||)^ - A - (fc 2 - P2)2][p| - 2t\qeB\ - to2] ’ 

(IIL6) 

and the numerator J\f^f {xi,yi', ^ 2 ,^ 2 ) defined by 
{Xi,yi]k2,P2) 

= k2,P2)(^e + mn)W^^'' (yi-,p2,k2)g^iv, 

(IIL7) 


Vhi«^(zi;4,£2) = (ni.8) 

To evaluate the numerator (xi,yi \ ^ 2 ,^ 2 ), we use the 

definition of Pn'^ from (III.261) . the properties of pro¬ 
jectors, = P^^ and P^'^P^'’ = 0, the relations 

from l|II.39|) and the following Dirac algebra 

7 II . 7_l = 0 , 

7 " • 7|| = 7-L • 7"^ = 2, 

7 '^“ 7 “ 7 mii =- 27 “^ 7 '^" 7 “ 7 Mr =-27“", (ni. 9 ) 

and arrive after some computation at 


{Xl,yi\k2,P2) 


= -2U 


+h 


iPi"^n,/+^^(a;i,fc2)/, ^%x,,p2)f, ^Hyi,P2)f+^^iyi,k2) 


+'P-^Tlnfj, ®‘*(a;i,fc2)//'’'’(a;i,p2)//'*'’(yi,P2)/„ *’(^ 1 ,^ 2 ) 


(cri, fc 2 )/ 7 '’(a:i,P 2 )/,+^’(yi,P 2 )/+^''(yi, fc 2 ) 


+iPi‘'^n„n,/+*‘'(xi,fc2)/r’(a^i,P2)/, ^\yi,P2)f-^^iyi,k2) 

'+^ fn"'’{xi,k2) {xi,P2)f^"'’ {yi,P2) {xi,P2)fi"'' {yi,P2)^ /,^'''’(yi, *2) 

+P^^''llnfn"'’ixi,k2) (llef["'’{xi,P2)f["’’iyi,P2) + f^"'’ixi,P2)f^"'’iyi,P2)^ fn"'’iyi,k2) 


r + S, 


(III.IO) 


Plugging at this stage {xi,yi;k 2 ,P 2 ) from (IIII.IOI) in 
(IIIL6I) . using P^'^ = P± and P^^^ = P^^ performing the 
integration over Xi and yi by making use of the method 
introduced in the Appendix, and evaluating the integrals 


using the standard Feynman integration method, we ar¬ 
rive after a lengthy but straightforward computation at 
[see also mn] 


h = 

Je = 


f d^PW _ I _ 

J (27r)3(p2-M2)((k||-p||)2-p'^)’ 


f d^PW _ 

J (2^)3 (p| 


g|| 

^|)((k|| - P ||)2 


P'i)’ 


(III.ll) 

I 


s£‘?)(A:) = + P-A^^^) + mq(P+d^^ + 


( 111 . 12 ) 
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with the coefficients given by 

OO 

4+) = -2eVE 

OO 

£^0 ' 
OO 

!)(+) = +2e^q^Y^ 
'<(+) 


<Pp'^ 1 

(27r)2 n!£!^ 
1 

(27r)2 n!£!^ 

1 

(27r)2 nlf! 

1 

(27r)2 n!£! 

V 

for positively charged particles, and 

!y^ f d^p'^_ 1 

£^7 (27r)2n!£! 

I f d'^p'j_ 1 

^7 (27r)2n!£! 


i=0 ' 

OO 

= +2e2g2^ 

OO 


4 +) = +2e‘^q^ 


aV-) 

= -2e^q- 


= -2e^q- 

£)(-) 

= +2e‘^q‘ 

Gj"^ 

= +2e^q- 

gI") 

= +2e‘^q‘ 


V [^±- 

^7 (27r)2n!£! 


f ^ 

£^7 (27r)2n!£! 

V [^±- 

^7 (27r)2n!£! 


for negatively charged particles. 


M, m) k||), 

“''nn„Gi(K; M', m')>7^(p'j_5 ky), 

-"ffl„n,G2(K;M",m",M"- l,m" - l)I,(pY,k||), 
-'‘[m,Gi(/t; M, m) + Gi(/c; M", m")]lKpl, ky), 
“'"n„[nffl£Gi(K;M" - l,m" - 1) + nGi(K; M', (p'j_) ky), 

■e“'"ffl„Gi(K; M, m) (p^, ky), 

■e“'^nn<!Gi(K; M', to') Jf(p'j_, ky), 

■e-«m„nfG 2 (K;M",TO",M" - 1,to" - l)If(p'^,ky), 

■e-''n„[£Gi(/c; M, to) + H^Gi( k; M", to")]I^(p'^, ky), 
■e“"[n£Gi(«;; M” — 1, to" — 1) + nn^Gi(K; M', to')]!^(p)l, ky ), 


(III.13) 


(III.14) 


In the above relations, ky = (fco, 0, 0, fca), p^ = kj^ — pj_ 

with kx = (0, fci,fc 2 , 0 ), Mj = 2(.\qeB\ +to^, k = 
with £q = \qeB\~^^‘^, and Gi(K; a,/3), z = 1,2 are defined 
by 

Gi(K;a,/3) = «:“"^[G 4 ^_^l(K)]^ 

G2(K;a,/3,7,<5) = k““'’[g4/3+i('«)]44i('^)]> 

(111.15) 


where G^(k) is the confluent hypergeometric function of 
the second kind [2^. Moreover, we have 


TO 

= min(n, £ — 1 ), 

M 

= max(n, £ — 

to' 

= min(n — 1 , £), 

M' 

= max(n — 1 , 

to" 

= T[mi{n,£), 

M" 

= max(n,t'). 


(III. 16) 

and 


2 ^^(Pj.:ky) 

^^(Pj.:ky) 


1 


1 


2 ^" /4k2M2 - {p'l - k| - M 2)2 

(p'i - ky - 


zarctan 


1 

47rk^ 


Y4kfM2-(p'i-k2-M2)2' 

P'i - k 2 + M 2 

Y4k|M2-(p'i-k|-M7)2' 


^4k2M2 - (p'i - kf - M 2)2 


arctan 


p'l - k 2 + M 2 


2 P'li 


(III.17) 


These integrals arise from Ii and Ji in (IIII.llI) . As afore- two-dimensional integration over p'j^ appearing in (IIII.13I) 
mentioned, the summation over Landau levels £ and the 





























and (IIII.14I) will be numerically performed in Sec. IV Al 
At this stage, let us only notice that in the LLL (n = 0) 
and from (IIII.13I) as well as 
and ^ from (IIII.14I) vanish. The general structure of 


for positively and negatively charged fermions in the 
lowest and higher Landau levels is therefore given by 


J 



0 , 

)4|| iV+A‘'f + iP_ aL®)) + mg{V+df + 


for n = 0 , q = + 1 , 

for n = 0, q = —1, (III. 18) 

for n ^ 0 , q = ± 1 . 


In Sec. El we will combine the above result from (IIII.18I) 
with the general structure of the free fermion propagator 
from (III.411) . and will derive the general structure of the 
one-loop corrected fermion propagator for nonvanishing 
magnetic fields at zero temperature. 


and determine ij and Jj from (IIII.19I) within this ap¬ 
proximation for = 0 [see Sec. IIII B 2| . 

1. First method: Exact results 


B. One-loop self-energy of magnetized fermions at 
finite temperature 


In this section, we will show that the general structure 
of for S ^ 0 and at T 7 ^ 0 is given by (IIII. 21) . We will, 
in particular, determine the coefficients A^£^ ^ 
and appearing in (IIII. 21) for i? 7 ^ 0 and at finite 
temperature T and chemical potential /i. 

To start, let us consider from (IIII. 61) with 

{^i:yF^ 2 ,P 2 ) presented in (IIII.10|) . To determine 
at finite temperature and chemical potential, all we 
have to do is to evaluate the integrals appearing in (IIII. lip 
at finite T and fj,. We thus focus only on the integrals 


-'t 


+ 00 

E 


n— — oo 


+00 

E 


n— — oo 


dp3 _ 1 _ 

27r (p2 - M2)((k|| - P||)2 - p'^) ’ 


dp3 _?||_ 

27r (p2 -M2)((k|| -P||)2 


P'i)’ 

(III.19) 


that arise from (IIII.Ill) with the standard replacement 


Using the standard Feynman parametrization, the in¬ 
tegrals ij and jJ from (IIII.19P are given by 

ij = i j dx j 

jJ = *y dx J + X 

(III.22) 


where Sg"'\uj) as well as are defined by 


+ 00 






'2\m 


n——oo 


(£ 2 -a; 2 )^’ 




n— — oo 


i 


(III.23) 


with f > 1 and m > 0. Moreover, = p^ + a;p'5_ — 2^(1 — 
a;)ky -I- (1 — x)Mg with Mf = 2^\qeB\ + m^. To proceed, 
we use 

= -^Nfico), (III.24) 


/ 



+00 

E 


n=—00 


/ 


dp3 

27r 


(III.20) 


with Nf{uj) = njr(w) + n^{uj) and Nf{uj) = n^(uj) — 
nj(uj), with the fermion distribution function 


Here, Pq = iujn — p with the fermionic Matsubara frequen¬ 
cies ujn = {2n+l)TrT. To compute these integrals, we will 
use two different methods. First, using the method pre¬ 
sented in Q, we will evaluate them without any approxi¬ 
mation [see Sec. IIIIB 1| . We will then perform the HTL 
approximation in a weak magnetic field, characterized by 
fca <C T as well as by 

< IpY \<^P3^T, (III.21) 


'(w) = 


1 


g/3(c^TM) + 1’ 

and the following relations for £ > 2 Q 


5f(u;) 


1 dS^°}^{uj) 

2{i — 1 ) 0 ; dio 

1 dsf\[pS) 
2 (£ — 1 ) 0 ; djjj 


(III.25) 


(III.26) 
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We arrive therefore at 


'52°^^) - + 


1 Nf{oj) N'f{uj) 


The general structure of Tin\k) at finite temperature and 
chemical potential is thus given by [see also (IIIL2I) ] 


(c.) = - 


4uj^ Auj^ 4uj^ 
4:0} 


(111.27) 


p_^(9)) +^ 3 ( 75 ^ 5 ^ + p_s(9)) + 1 ^ 3 _d(‘?) + TO,(iP+C7^«) + (111.28) 


with the coefficients given by 


yl^+) = -2eV£ J «mfGi(/t;M,m)77/°^^(pl,k||), 


= -2e2'-2 


e=o • 

°° - j 2 „/ 


”nn„Gi(K;M',m')j/°^^(p'j.,k||), 


5^+^ = - 26 ^ 


e=o • 


9 "^/ ”m^Gi(At;M,m)J-/"^^(pY,k||), 


£=0 ' 

^ - 72 ^/ 


= -2e2g2y^ [ ''nn„Gi(K; M', m')77/^^^(p'j_, ky), 


^ J ( 27 r)^ nh?! 

^=0 ^ 

D<«.+2»vy'/54 ‘ 


V2.^ / 72^-j7ie-Vn„n,G2(K;M",m",Ar-l,m"-l)lJ'(p'j.,k|| 

I TTj ■ 


cf) = + 2 e 2 g 2 ' 


/* r /2 / 1 

.+ _ , = ^ ^ e-«[ffl,Gi(Ac;M,m) + Gi(«;M",m")]lJ(pl,k||), 

ci+) = +2eV^ / -^—e-«n„[nm,Gi(K;M"-l,m"-l) + nGi(At;M',m')]Xj(pl,k 
i=o I ^1 ^- ■ 

for the positively charged particles, and 

= -26^9^^ / 4%^e"'"m„Gi(K;M,TO)j/°^^(p(L,k||), 

= -2e2(?^y^ / 4^T^e"'‘nn£Gi(«;;M',m')J'/°^^(p(L,k||), 


11 ), (III.29) 


“7 (27r)^ n!£! 


= -2e^q^ 


^ '‘m„Gi(/c;M,m)J'/^^^(p(L,k||), 


fco 


sy* = - 2 evf 


e-“nn,Gi(K;A/',m')jf’^(pY,k||), 


J (27r)^ n!^! 

£=o ^ ' 


g(-) = +2eVE / ^-|^e-"m„n,G2(K;;M",m",M"-l,m"-l)lf(pY,k||), 

00 ^ 1 

Gt^ = +2eVE J ^;^e-"n„[^Gi(K;M,m) + n,Gi(«;M",m")]lJ(pi,k||), 


Ci"^ = +2e2g2'^ / ^^-Le-"KGi(K;M"-l,m"-l) + nnfGi(K;M',m')]lJ(pl,k||), (III.30) 


(27r)2 n!^! 
^=0 ^ 


for negatively charged particles. 


In the above relations, k is again given by k = ^ , aswellasM, M and M are presented in piI.16[ ). More 

Gi{K;a, f3),i = 1,2 are defined in (IIII.15I1 . and m, m', m" 
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over, lJ(p'j_,k||) and ky),z = 0,3 are given 

by 


^{0)T 

'Jl 

7-(3)T 

'Jl 


Jq J ^ 


dx 


dps 

27r 




ko 


j ^ ^ (III.31) 


with 52 °^(w) and from (IIII.27I) . Let us notice 

that, for fj, = 0, because of , we obtain 


for both g = +1 and g = -1. As in T = 0 
case, the summation over Landau levels and the two- 
dimensional over p^, appearing in (IIIL29I) and (IIIL30|) . 
as well as the integration over and x, appearing in 
(IIIL31|) . are to be numerically performed. Here, similar 
to the T = 0 case, in the LLL (n = 0), all the coeffi¬ 
cients A^\B^\D^~'> and for negatively charged 
particles vanish. Similarly, for positively charged parti¬ 
cles, the coefficients and vanish for 

n = 0. The general structure of in the lowest as well 
as higher Landau levels is therefore given by 


J 



V+ )) , 

MB+A^+ + V-A^^^) +^ 3 (^+ 5 +^ + r-B'f'>) 


for n = 0 , q = - 1 - 1 , 
for n = 0, q = — 1 , 
for n 7 ^ 0 , q = ± 1 . 

(111.32) 


In Sec. HYi we will combine these results with the gen¬ 
eral structure of the free propagator from (111.411) . and 
will derive the general structure of the one-loop corrected 
fermion propagator for nonvanishing magnetic field at fi¬ 
nite temperature. 


2. Second method: HTL approximation 

In this section, we will evaluate the integrals (IIIL19I) . 
using the modified HTL approximation (IHL21I1 . For sim¬ 
plicity, we consider the /i = 0 case. Let us introduce, as 
in Sec. Ill A1 

Ef ^ pI+MI 

El = (p3-fc3)2 + p'^, (IIL33) 

and rewrite the integrals appearing in (IHL19I) as 

^J(p'j.:k||) = * J ^ 2f^b3;p'j_,k||), 
JUp'x^'^w) = i J ^[7oS^°kP3;p'_L,k||) 

-73P3E;f^(p3;p'_L,k||)], (111.34) 

with 

+00 

5^°^ - ^ X! *w„A/(za;„,£;i)At,(fco - za;„,£’2), 

n— — oo 
+00 

^ X! ^fi'>‘^ri,Ei)Ab{ko - iuJn,E2), (111.35) 

n— — oo 

similar to (IH.5|) . with Af{po,Ei) and At,(/co — po, ^^ 2 ) de¬ 
fined in (|H.3I) . and Ei,i = 1,2 given in (IHI.33I) . The 
summation over Matsubara frequencies can be performed 


using the same method as described in Sec. Ill Al Similar 
to (IH. 6 I) . we therefore have 

7 ^(0) _ 1 + /b(s2L'2) — //(siL^l) 

^ , 4£/2 fco — siEi — S 2 E 2 

Sl,S2=± 

.^(3) _ _ 5iS 2 1 +/&(52A2) — //(siAi) 

^ , 4:EiE 2 ko — siEi — S 2 E 2 

Si,S2=± 

(111.36) 

with the distribution functions ff/b{E) defined in (|H.7I) . 
and Ei, i = 1,2 given in (jHL33p . Under the modified HTL 
approximation introduced in (IHL21I) . we set 


Al 


Eo 


P 3 


Ml 

2P3’ 


P 3 - ^3 + 


/2 

P± 

2p3 ' 


and have consequently 


(111.37) 


//(Al) 

ME 2 ) 


flips) 


Ml dffjpo) 

2p3 dp3 


fbips) + 




dfbjps) 

dps 


(111.38) 


Using furthermore the relation //(—w) = 1 — //(w) and 
fb{—ui) = —1 — fbiuj), we arrive after some computation 
at 
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with 


jT ^ 

Jj = 


(111.39) 














































87rT2l 0 


+ 


rT(2) I t(2)i 2fc|-(2) M^2 

^ih + 1 / J 


2k3T 


(fcg + /c|)(p'i-M^2) ■ 

-. 2\2 ^ 


(k^ 


f(4) 


r(fc2 + A:|)(p'l-Mg) 1 

l(3) 

\ikl + kl){p'l-M^) , J 

[ (kf)^ i 


[ (kp^ J 


r(3) 


+ 




/2 
P J. 

2 T 2 
1 


(fcg + fc|)(pl-M2) 


1 


f(4) , *3 

■■■fc “T 


-2fc2p'l + M2(fc2 ^ fe2) 


- 1 


J(3) 


■fif 


-If] 


(k^)= 

_ (2) , Mfhip'l - M!) =is) 

2Tkf,''f 2T2(k|)2 f 




jf^(p\,k||) = 


1 r T’rid) , i(i)i , fc3(p'l-M2) 

AnTl k|^^^ +^/ (k2)2 ^ 

r p'l(fcg + fcg) -2fcgM^2 -| fc3pj(p'i -Mg) -(3) 

2T(k2)2 J b 2T2(k2)2 b 

^g+fcg)(p'i-M 2 ) 


+ ^I( 1 ) _ 
kg f- 


-If] 


8nT \ fc 3 k| ** 

, 1 r(fcg+fcg)(p'i-M 2 ) ^ 

“T T /I o\2 


2 T2(k2)2 

M?-(^2) ^|( 1 ) 

' f ^ kg ^ 2T2fc3 


Tkg"/ 


(kS) 


- 1 


t(3) 


+ 


/2 
P J. 

2 r 2 fa 


(k|) 

(fcg + fcg)(p'g_-Mg) , ^ 

ru2\2 ^ 


l( 2 ) , 1 

^ +fc3 


(fcg+ A:g)(p'g^-Mg) ^ g 

(kg )2 


r(2) 


f 


-2A:gp'g_ + Mg(fcg + fcg) 


In the above relations 


r(") 


(-) = r^- 

y 

lif W = / $/b//(2/)> 

*/2 y 

!!;}(-) - 


dy dfb/f{y) 

2/” dy 


with z = , and 


■^b( 2 /) gy _ 2 ^ ’ //(d) gy _|_ 2 ■ 


(111.41) 


(111.42) 


The general structure of (fc) for hot and magnetized 
fermions in the above HTL approximation is therefore 
given by (IIII.32I) with the coefficients and 

from (|III.29I) as well as l|III.30p . and the integrals Xj 
as well as / = 1,2 given in (IIII.40I) . 


as a special case of T ^ 0 case, because by comparing 
from (|III.18I) for T = 0 with (|III.32I) for T 7 ^ 0, it turns 

out that for T = 0 can be determined from (IIII.32I) 
by setting = vlf, i.e. by removing the anisotropy in 
the components o fk|| = (fcp, 0,0, ks), appearing in (IIII.32I) 
in comparison to (IIII.18I) for n = 0 and n 7 ^ 0. In this sec¬ 
tion, we are not interested on the numerical values of the 
coefficients. Our main goal is to use the general structure 
of the dressed propagator of hot and magnetized fermions, 
and, performing an analysis similar to what is presented 
in Sec. (|II Al) . to determine the properties of possible ex¬ 
citations arising from the poles of this propagator. 

To start, let us therefore consider the series expansion 
of the full fermion propagator in the momentum space 

sL^\k), 


(/) = (/) + (/)s(g)5(«) + • •. . (IV.l) 


IV. GENERAL STRUCTURE OF THE DRESSED 
PROPAGATOR OF HOT AND MAGNETIZED 
FERMIONS 

In (|II.4ip . the free propagator of magnetized fermions 
is presented in the momentum space. The general struc¬ 
ture of the one-loop self-energy of these fermions at finite 
temperature is presented in (jIII.32ll . In this section, we 
will combine these two results, and determine the one-loop 
corrected propagator of hot and magnetized fermions up 
to one-loop level. The case T = 0 will be then considered 


Truncating this series after the one-loop contribution in 
the second term on the r.h.s., and using (III.41|) as well 
as (IIII.32I) . the general structure of Sn ^ of magnetized 
fermions at T 7 ^ 0 in the lowest (n = 0 ) and higher {n 7 ^ 0 ) 
Landau levels up to one-loop level reads 



■p 


¥ 11 -^ 






for n = 0, q = - 1 - 1 , 

for n = 0, q = —1, (IV.2) 

for n ^ 0, q = ± 1 , 
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where 


with the numerator' 




= l^o(V+A^^^ +V-A^^^) +iP_Bi‘'^) 

+ mg{V+df + (IV.3) 


= hil-B^^^), 

= 771,(1 + 

= s,v'2n|geS|(l-£»(«)), (IV.5) 

Sn'^ for n = 0 and n ^ 0 can be simplified as"^ 

^+(^70 0 ^. -73 0 + +C+ J 


\{i)\ 


(9) 


5, 


with 


■jyT+V 

7’^(7ofco—73fc3+mg) 


•D, 


rr' 


for q = +1, 
for g = — 1, 


(IV.6) 


p(+) = 


(+)2 . (+)2 , (+)2 


and 


5 ( 9 ) (-fo _ (^) 


’ 


(IV.7) 


(IV.8) 


^ The 7 -matrices in the Dirac representation are given by 
10\ / 0 cr\ f 0 1 


70 = 


0 -1 


7 = 


-<T 0 


75 = 


1 0 


->(9) 

-'3 


= 7o(c(«V++4«V_)+72^ 
+73(c1'V+ + 4«V_) 

+ (c^'^So 72 + 7273 foCs'^^ 707273) S3 


are obtained from (jIII.32ll . At T = 0, Sn'^ is given by 
(IIV.2P and (IIV.3P with for g = ±1. The poles 

of the fermion propagator (IIV.2I) can be determined by 
computing 

det(Jt|| — rUg — Sq’*’^), for n = 0, q = -1-1, 

0 = < det()t|| — mg), for n = 0, q = —1, 

det(^„ - mg - T,n'^), for 77 ^ 0, q = ±1. 

(1V.4) 


+ci’^'>V++c[l^V-, (IV.9) 

and the denominator 

= ^(9)4 ^ 2d(«)2 

+ . 

(IV.IO) 

In (|IV.9I) . the coefficients f = 1, • • • ,10 are given by 


This will be done numerically in Sec. IV Al for T ^ 0 [hot 
magnetized QED plasma], and in Sec. IVBl for T = 0 
[cold magnetized QED plasma]. To study the proper¬ 
ties of possible excitations arising from the pole of the 
dressed fermion propagator, let us consider Sn'^ from 
(IIV.2p . Defining 


ry<J) — 
'-3 — 




- - 

,( 9 ) 2 ) 


5 

( 

[ a'f 

_ 6 ( 9)2 _ 

C^) 


5 


2 + cf«)c( 9 ) + 6 

,(9)5(9) 



i,“( 

'5(9)2 


,( 9 ) 2 ) 

+ 6 L^ 7 ( 9 ) 2 , 


i,7( 

'5(9)2 

- a ^ f - 

cV^)^) 



(c<-> 






(c? 

5(9) _ 





(<■? 

5(9) _ 





c“( 


- - 

,( 9 ) 2 ) 



c7( 


_ 6 ( 9 ) 2 . 



(IV.ll) 


f'- 9 ; _ 

'-'10 — 


In what follows, we will determine the eigenvalues and 
eigenvectors of the numerators of Sn'^ for two special 
cases: 

Case 1: n = 0, q = ±1, mg yf 0, 

Case 2: 0, q = ±1, mg = 0. 


Case 1: Properties of massive and positively charged 
fermionic excitations in LLL 

Let us consider the one-loop corrected propagator of 
positively charged massive fermions in the LLL from 

(ITvel) 


s^-^\k) = 


iP+(7oaV^^-736V+^+cV+^) 


n 


(+) 


, (1V.12) 


and S = diag(cr.<T). 


Here, the relations + V— = 1 and P-f — P_ = E 3 are used. 
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with the coefficients given in (IIV.5|) and the denominator 
given in (IIV.7I) . The eigenvectors of the numerator 
of are 


4+^ 

44 


( 0 , 0 , 0 , 1 ), 
( 0 , 1 , 0 , 0 ), 


1 




+ 


jW2_5(+)2^0,&^+\o^ 


(IV.13) 


with the eigenvalues 


The eigenvectors from (IIV.17I) satisfy 


r+44 

= 44^ 

or 

^^44 

= +44 

v+4~^ 

= 42~\ 

or 


= +' 02 ”^ 

v-4 ^ 

= 44^ 

or 


= -03”^ 

r-44 

= 44^ 

or 

^344 

= -04”^ 


(IV. 19) 


Similar to the previous case, only 4 = 3,4 corre¬ 

sponding to negatively charged fermions with negative 
spins (spin down) in the LLL, are acceptable. In con¬ 
trast to the case discussed in Sec. Ill A1 44^^ = 3,4 

and 4 = 3,4 are neither eigenvalues of the helicity 

nor those of the chirality operators. These operators are 
defined in (III.171) and pi.181) . respectively. 


|o,0,c^+) 

- 4^4^ 


.c4 + ^ 

44^ - 

b4"\- 






(IV. 14) 

The eigenvectors from 

(IIVT^ 

satisfy 



v-44 

= 44. 

or 

^344 

= -44. 

1 

v-44 

v+44 

= 44. 

= 44. 

or 

or 

^344 

'^344 

= -44. 
= +44. 

’ (IV.15) 

1 

v+44 

= 44. 

or 

'^344 

= +44- 


From the above four 

eigenvectors only 44'^. 

1 = 3,4 


are acceptable. They correspond to positively charged 
fermions with positive spins (spin up) in the LLL. 


Case 1: Properties of massive and negatively charged 
fermionic excitations in LLL 

According to (IIV. 6 I) . the one-loop corrected propaga¬ 
tor of negatively charged and massive fermions does not 
receive any contribution from the one-loop fermion self¬ 
energy 


= (IV. 16 ) 

T)^ J 

Here, 1^4 is given in (IIV.7|) . The eigenvectors of the 
numerator of 5g ^ are 

44 = ( 0 , 0 , 1 , 0 ), 

44 = ( 1 , 0 , 0 , 0 ), 


4^ '’ = 4 ( 0 , |k||| - fco,o, 1 ), 

^3 


44 = -r (o,-mil “ ^ 0 , 0 , 1 ), 

(IV.17) 

with the eigenvalues 


{ 0 , 0 ,m, - |k|||,TO, -f |k|||}. 

(IV.18) 


Case 2: Properties of massless fermionic excitations in HLL 


The one-loop corrected fermion self-energy of mass¬ 
less fermions for n yl 0 and q = ±1 is given by (I1V.8I) - 

PV.III) with 4^ = 0. To study the properties of possible 

fermionic excitations, let us simplify s4 , in analogy to 
the results presented in Sec. Ill A1 as 


s44) 


444) ^ 44 Ck) 
2v4 4) 2v4 4) ’ 


(IV.20) 


with the numerators 

44 = MV+ia4 - b4) + V-{a4 + b^)] 

-2d^‘^h2VR + iz[V+{a4 - b4) - V-{a4 + 6^^)], 
M4 = lc^[V+{a4 + b4) + V.{a4 - b4)] 

-2S^h2VL - 4V+{a4 + b4) - V-44 - b4)], 

(IV.21) 

and the denominators 

V^/ii = {44 ± b4){a4 T b4) - (IV.22) 

Nontrivial eigenvectors of the numerators J^4 -^4 

are given by 

Vf = (1,0,0,0), 

= (0,1,0,0), (IV.23) 

and 

Vf = (0,0,1,0), 

= (0,0,0,1), (IV.24) 

respectively. They have trivial eigenvalues, and satisfy 


VL+t 

= 4. 

and 

SaVf = 

+V[, 

VlV^ 

= Vi, 

and 

^3Vi = 

-Vi, 

VrV^ 

= Vf, 

and 

SaVf = 

+Vi 

-PR+i 

= Vi, 

and 

V.3Vi = 

-Vi 
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As yet the above results (IIV.23I) - (IIV.25|) are independent 
of the choice of electric charges q. They are valid for both 
positively and negatively charged particles, and can thus 
be summarized as 


,(9) . J 

r vf 

with 

q — +1, 

S3 

= ±1, 

X = 

-1, 

^ -i 

l vf 

with 

q = -1, 

S3 

= ±1, 

X = 

-1, 

1(9). J 

^ vf 

with 

q = +1, 

S 3 

= ±1, 

X = 

+1) 

* -j 


with 

q = -1, 

S3 

= ±1, 

X = 

+1. 


(IV.26) 

Here, S 3 and y denote the eigenvalues of the spin and chi¬ 
rality operators, S 3 and Vl/Rj respectively. Let us notice 
that the appearance of two energy branches, arising from 
the poles of two denominators and 2?^^ in (IIV.20I) . 
can be regarded as an evidence of the appearance of ad¬ 
ditional fermionic excitations. In the next section, we 
will study the spectrum of Dirac particles at finite tem¬ 
perature and for nonvanishing magnetic fields. We will 
show, that in the limit of soft momenta and weak mag¬ 
netic fields, new excitations appear, which will be referred 
to as hot and magnetized plasminos, in analogy to the ex¬ 
citations appearing at finite T and vanishing B. 

V. NUMERICAL RESULTS 

In this section, we will numerically solve the energy 
dispersion relations arising from the one-loop corrected 
fermion propagator, Sn \ for a number of special cases. 
In Sec. IV A[ we will first consider at nonzero T and 
B in the lowest (n = 0) and higher (n ^ 0) Landau levels, 
separately. We will focus on both massive and massless 
case. In Sec. IV Bl we will then determine the spectrum 
of massive fermions in a cold and magnetized electromag¬ 
netic plasma. 

A. Plasminos in a hot and magnetized QED plasma 

In the previous section, we have analytically deter¬ 
mined the general structure of the one-loop corrected 
fermion propagator in lowest and higher Landau lev¬ 
els. For n = 0, Sn^ is given in (IIV. 6 |) . and for n ^ 0 in 
(IIV.8|) - (|IV.11I) as well as in (IIV.20|) - (|IV.22I) for the special 
case of massless fermions. Except for negatively charged 
fermions in the LLL, the one-loop corrected fermion prop¬ 
agator is, in particular, given in terms of nontrivial coeffi¬ 
cients and , which are defined in (I IV. 51) , 

and are given in terms of the coefficients , B ^£^, 
and , arising from the one-loop self-energy of charged 
fermions. In Sec. IIII Bl we have analytically determined 
\ and at finite T and H, using two 

different methods: In Sec. im B II the exact expres¬ 
sions of these coefficients at finite T and for nonvanish¬ 
ing B are presented in (IIII.291) and (IIII.301) , with if and 

J = 0, 3 from (IIIL3I|) . and in Sec. IIII B 21 they are 


evaluated using a HTL approximation in a weak magnetic 
field, and are given by the same (IIII. 291) and (IIII.301) . with 
if and f = 0, 3 from (IIII. 401) . 

In order to explore the spectrum of fermionic exci¬ 
tations from the poles of Sn'^ in the lowest and higher 
Landau levels, the coefficients and 

are first to be determined as functions of kg and k^. To 
do this, we have numerically evaluated the integration 
over p'j^, x andps appearing in (IIIL29I) - (IIIL31I) . as well as 
the integration over y appearing in Iq , I^// and Ib/f from 
(IIII. 411) for a large number of fixed kg and k^. The sum¬ 
mation over Landau levels appearing in (IIIL29I) - (IIIL30|) . 
has also been performed numerically. In this way, it was 
possible to find the best fits for and 

as functions of ko and k^. We have then considered the 
fermion propagators in LLL (n = 0) and HLL (n = 1) sep¬ 
arately,® and solved the energy dispersion relations, aris¬ 
ing from the poles of the propagators. In what follows, we 
will only report the corresponding numerical results for a 
number of special cases. 


1 . Special case: T 7^ 0, n = 0 for rriq 7^ 0 and rriq = 0 

Let us start with the dispersion relation of a massive 
and positively charged fermion in the LLL. It can be ob¬ 
tained by numerically solving det()t|| — rriq — = 0 

from (I IV. 41) or ik) = 0 with from (IIV.7I) . In 
both cases, it reads 

kl{l - - ^,(1 + = 0 . 

(V.l) 

As concerns the negatively charged massive fermions, the 
corresponding dispersion relation arises from (Dq ^ = 0 , 

kl = kl + mf (V.2) 

Here, the definition of from (IIV.7I) is used. The 
numerical results are demonstrated in Fig. [31 In Fig. 
|3]Ja), the k'i/T dependence of k^/T is plotted for positively 
(thick blue curve) and negatively (thin red curve) charged 
massive fermions in the LLL (n = 0) for b = eB/T^ = 5, 
^ = rriq/T = 0.25 and in the regime k^/T > 0.2.^® To 

do this, we have used (|HL29I) with if and i = 0,3 

from (|HL31I) . and numerically determined the dependence 
of A^\ and appearing in (IV.II) . on kg and kg. 

This results in the thick blue curve in Fig. |3Ka). The 
thin red curve in Fig. |3Ka) is the positive energy branch 
of the dispersion relation of negatively charged massive 


® For HLL, we considered only n = 1 case. All n > 1 cases can be 
evaluated in a similar way. 

At finite temperature, the free parameter eB and mq are scaled 
with T. 
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0.00 0.02 0.04 0.06 0.08 


k3/T 

FIG. 3: (color online), a) The exact ks/T dependence of ko/T 
for positively (thick blue cnrve) and negatively (red thin cnrve) 
charged fermions in the LLL (n = 0) for b = eB = 5 and 
^ = mqjT = 0.25, in the regime fes/T > 0.2, which arises 
from (IVll) for q = +1 and (T^ for q = -1. b) The HTL- 
approximated ks/T dependence of ko/T for positively (thick 
blue curve) and negatively (thin red curve) charged fermions 
in the LLL (n = 0) for b = ^ — 0.5 and ^ ^ = 0.25, in 

the regime ks/T < 0.2. 


fermions from (IV.2I) for 6 = 5 and ^ = 0.25.^^ The regime 
ks/T < 0.2 is then explored using and 

from (IIIL29|) with if and ,i = 0, 3 from (|IIL40I) for 
b = 0.5 and ^ = 0.25. In this regime the HTL approxi¬ 
mation described in Sec. IIII Bl might be reliable. Similar 
to Fig. EKa), the thick blue and thin red curves corre¬ 
spond to the positive energy branches of the dispersion 
relations (IV.II) and (IV.2|) . respectively. Let us notice at 
this stage, that according to the results from (IIV.15I) and 
(IIV.19I) . and the subsequent descriptions, the blue and 
red curves in Fig. [3] are the energy dispersion relations of 
positively and negatively charged massive fermions with 
positive and negative spins. 


The negative energy branches of the energy dispersion relations 
are not plotted in the figures demonstrated in this section. 


a) T:/:0, n=0, b= 0.5, ^=0 



ks/T 



FIG. 4: a) The exact ko/T dependence of ko/T for positively 
(blue curve) and negatively (red squares) charged fermions in 
the LLL (n = 0) for b = eB/T^ = 0.5 and ^ = mq/T = 0, 
in the regime ko/T > 0.2, which arises from (El for q = -1-1 
and (IV.4II for q = — 1. b) The HTL-approximated ko/T de¬ 
pendence of ko/T for positively (thick blue curve) and nega¬ 
tively (thin red curve) charged fermions in the LLL (n = 0) 
for h — eB/T^ = 0.5 and ^ = mq/T — 0, in the regime 
ko/T < 0.2. 

Comparing the energies of fermions with positive and 
negative spins from Fig. EJa) [exact results] with the 
corresponding energies of these fermions from Fig. Hb) 
[HTL-approximated results], it turns out that in the 
regime k^/T > 0.2, their T-scaled energies, kq/T, increase 
with increasing k^/T, while in the regime k^/T < 0.2, the 
energy of fermions with positive spins (thick blue curve) 
decreases with increasing k^/T, in contrast to the energy 
of fermions with negative spins (thin red curve). 

Performing the above analysis for the case of massless 
fermions in the LLL, it turns out that the small difference 
between the energies of fermions with positive and nega¬ 
tive spins, appearing in Fig. da) disappears in the limit 
of vanishing fermionic mass. This can be observed in Fig. 
da), where the exact k^/T dependence of kq/T is plotted 
for n = 0, 6 = 0.5 and ^ = 0 in the regime k^/T > 0.2. 
Here, the energy dispersion relations 

eo{l-Alf-kl{l-BXf = 0, (V.3) 


a) 7:^0, n=0, b=5, ^=0.25 
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kzfT 

FIG. 5: (color online). The exact k^/T dependence of fco/T for 
positively (panel a) and negatively (panel b) charged fermions 
for n = 1,6 = eB = 5 and ^ = mqjT = 0 in the regime 
fcs/T > 0.2. The results arise from the dispersion relations 
(IV31l and (IV.61) for left- and right-handed massless fermions. 
Blue curves and red squares denote the solutions of = 0 
and = 0, respectively. 



0.0 0.2 0.4 0.6 0.8 1.0 

/cs/T 



0.0 0.2 0.4 0.6 0.8 1.0 

/cs/T 


FIG. 6: (color online). The exact fcs/T dependence of feo/T 
for left-handed (panel a) and right-handed (panel b) fermions. 
Thick blue and thin red curves denote the solutions of = 

0 and 2?^^^ = 0, for positively and negatively charged particles 
q = ±1 for n = 1, 6 = eB/T^ = 5 and ^ = ra^jT = 0 in the 
regime fcs/T >0.2. 


for positively charged fermions, and 

kl = kl (V.4) 

for negatively charged fermions, are used. In the regime 
k^/T < 0.2, however, the qualitative behavior of the k^/T 
dependence oikq/T for massless fermions is similar to the 
corresponding results for the massive fermions [see the 
plots in Fig. HKb) and compare them with the plots in 
Fig. |3Kb)]. Let us notice, that although the fermionic ex¬ 
citations, whose energy dispersion relations are plotted in 
Figs. Hb), are massless, their rest masses (for k^/T = 0) 
are nonzero [see, in particular, the rest mass of the pos¬ 
itively charged massless fermion from Fig. fflb)]. This 
effect is related to the magnetic catalysis which is 
characterized by a dynamical generation of mass in the 
presence of very strong magnetic field, where the system 
is dominated by the LLL. As it is well-known, the dynam¬ 
ical mass arises from perturbative corrections to fermion 
propagator. According to the same argument, the fact 


that the negatively charged massless fermions have a zero 
rest mass [see the thin red curve in Fig. HKb)] is therefore 
related to the fact that these particles do not receive any 
contribution from the one-loop self-energy in the LLL. 


2. Special case: T ^ 0, n = 1 and mq = 0 

The one-loop corrected fermion self-energy of massless 
fermions in HLL is analytically computed in the previous 
section [see (jIV.20|) - (IIV.22l) ]. In this case, in analogy to 
the case of vanishing magnetic fields, described in Sec. 
Ill A1 two denominators appear for each g = -|-I and 

g = — I in (IIV.20I) . We therefore expect two different en¬ 
ergy branches for each g. The energy branch arising from 
= 0 corresponds to a left-handed positively or nega¬ 
tively charged particles, possessing both spin orientations. 
Similarly, the energy branch arising from 2?^^ = 0 corre¬ 
sponds to a right-handed positively or negatively charged 
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a) Ji^O, q=+1, n=1, b=0.5, ^=0 



/C3/T 



ks/T 


FIG. 7: (color online). The HTL-approximated fca/T depen¬ 
dence of ko/T for positively (panel a) and negatively (panel b) 
charged fermions for n = 1, 6 = eB/T^ = 0.5 and ^ = ra^/T = 
0 in the regime fcs/T < 0.2. The results arise from the dis¬ 
persion relations (IV.51) and (IV.61) for left- and right-handed 
massless fermions. Thick blue and thin red curves denote the 
solutions of = 0 and = 0, respectively. 



k^/T 



FIG. 8: (color online). The HTL-approximated ks/T depen¬ 
dence of ko/T for positively (panel a) and negatively (panel b) 
charged fermions for n — l,b = eB— 0.5 and ^ = niq/T = 
0 in the regime fca/T > 0.2. The results arise from the dis¬ 
persion relations (IV31l and (IV. 61) for left- and right-handed 
massless fermions. Thick blue curves and red squares denote 
the solutions of = 0 and = 0, respectively. 


particles with positive and negative spins. The explicit 
expressions for the energy dispersion relations of these 
particles are given by 

[fco(l-AL"^) + fc3(l-Bi''))] 

X [fco(l - - ks{l - - 2n\qeB\{l - 

= 0, (V.5) 

and 

X [fco(l - 71^^) + fc3(l - B^^'>)] - 2n|geB|(l - 
= 0. (V.6) 

They arise from = 0 and 2?^^ = 0, with 
defined in (IIV.22I1 . To solve these relations, let us first 
consider the coefficients A‘'^\ B^'^ and from (1111.291) 
for g = -1-1 and (IIII.30I) for g = — 1 with if as well as 
= 0,3 from (IIII.31I) . This yields the exact ks/T 


dependence of ko/T for massless left- and right-handed 
fermions. For n = l,b = eB/T^ = 5,^ = rriq/T = 0, 
the results are demonstrated in Figs. EDa) andEKb), in 
the regime k^/T >0.2. They correspond to g = -1-1 and 
g = — 1, respectively. As it turns out, there is no differ¬ 
ence between the solutions of = 0 (blue curve) and 
= 0 (red square) for positively charged particles. 
Similarly, the k^/T dependence of ko/T for negatively 
charged left- (blue curve) and right-handed (red-square) 
massless fermions are identical. 

To compare the energy dispersion relations of massive 
fermions for g = -|-1 and g = —1, we have plotted in 
Fig. [6] the exact k^/T dependence of ko/T for massless 
left- (panel a) and right-handed (panel b) fermions. Thick 

blue and thin red curves denote the solutions of = 0 
and = 0, for n = 1, 6 = 5 and ^ = 0 in the regime 

ko/T > 0.2. As it turns out, the energy increases with 
increasing ko/T > 0.2. Moreover, the energy of nega- 
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tively charged left- and right-handed massless fermions 
are smaller that their positively charged counterparts. 

The regime fca/T < 0.2 is explored in Fig. [3 where 
the HTL-approximated solutions of the energy disper¬ 
sion relations (jV.5l) for left-handed and (jV.Op for right- 
handed massless fermions are demonstrated. To deter¬ 
mine these solutions, the coefficients and 

from (IIII.29I) and (IIII.30I) with the HTL-approximated if 
as well as i = 0, 3 from (IIII.40I) are used. They are 

first numerically computed as functions of fco and ks for 
n = 1,6 = eB/T"^ = 0.5 and ^ = 0. Plugging then the 
resulting expressions in (IV.51) and (IV.61) , these energy dis¬ 
persion relations are numerically solved. The thick blue 
and thin red curves in Figs. [Tja) and [Tfb) denote the 
solutions of Vf = 0 and = 0, respectively. As it 
turns out, in the regime < 0.08T, the energy branches 
for the left- and right-handed positively [Fig. [3a)] and 
negatively [Fig. I3b)] charged massless fermions are split 
up, in contrast to the results for fca > 0.08T from Figs. 
|7| and [51 where it is shown that left- and right-handed 
massless fermions have the same energies. In Fig. jS] we 
have used the above HTL approximation to determine 
the k^/T dependence of ko/T for left- and right-handed 
massless fermions in the regime k^/T > 0.2. Here, we 
have chosen the same free parameters, n = 1, 6 = 0.5 and 
^ = 0, as in Fig. 13 The blue curves and red squares in 
Fig. 13 denote the solutions of — 0 and = 0, re¬ 
spectively. The results from Fig. |3 coincide qualitatively 
with the results from Fig. jS] where the exact solutions of 
(IV.5p and (IV.6p are demonstrated for n = 1,6 = 5 and 
^ = 0 in the same k^/T > 0.2 regime. 

Let us notice at this stage, that the appearance of 
two different energy branches in Fig. |3is in great resem¬ 
blance with the appearance of two energy branches for 
fermionic particle and plasmino excitations at finite tem¬ 
perature and zero magnetic fields. The latter case is dis¬ 
cussed in Sec. Ill A1 where it was shown that two energy 
branches in Fig. |3 correspond to fermionic excitations 
with positive (particle) and negative (plasmino) helicity 
to chirality ratios. At finite temperature and for nonzero 
magnetic fields, discussed in the present section, however, 
the massless fermionic excitations, whose energy disper¬ 
sions are demonstrated in Fig. |3 are, in contrast to the 
ordinary particle and plasmino modes at finite T, eigen¬ 
states of the chirality and the spin operators, T^l/h and 
S3 [see Sec. lYI for a proof]. Despite this difference, and 
because of the similarity between the production mecha¬ 
nism of new excitations in the case of nonvanishing T and 
B with the mechanism leading to particle and plasmino 
excitations at finite T and zero B, the new excitations will 
be referred to as hot magnetized plasminos. The crucial 
point is that, according to our findings, they seem to ap¬ 
pear only in the limit of soft momenta k^ T and weak 
magnetic fields eB <C where a HTL approximation is 
reliable, and where apart from LLL, higher Landau levels 
are to be taken into account. 

In the rest of this section, we will show that similar 
excitations appear at zero temperature and for nonvan¬ 


ishing moderate magnetic fields. 


B. Plasminos in a cold and magnetized QED plasma 

As we have described in the previous section, accord¬ 
ing to the notations used in the present paper, the zero 
temperature case can be regarded as a special case of 
the finite temperature case. This is because the similar¬ 
ity/difference between the general structure of one-loop 
self-energy of fermions at zero and nonzero temperatures 
from (IHL18I) and (IHL32I) . Thus, in the zero temperature 
case, the one-loop corrected propagator of fermions in the 
presence of a constant magnetic field is, as before, given 
by (IIV.6P - (|IV.7I) for n = 0 and (IIV.8I) - (IIV.11I) for n > 1. 
The only difference with the finite temperature case is 
that the coefficients and appearing in 

these equations, are to be redefined as 

^ ko{l-A^£), 

^ k^{l-A^f), 

= 7nq{l + &^^), 

= s,v'2n|geH|(l-£>(«)), (V.7) 

with and presented in (IHL13|) for g = -hi 

and (IHL14I) for q = —1. Similar to what is performed 
in the previous section, we have determined numerically 
the latter coefficients for a large number of /cq and k^. In 
this way, we were able to find the best fits for 
and as functions of fco and fca. Plugging then the 
corresponding expressions in a ±.'^, 6^^, and from 

(|V.7I) . we have numerically solved the energy dispersion 
relations for cold fermions in a magnetized QED plasma. 
In what follows, we will separately consider two different 
cases of n = 0 and n = 1 for massive fermions, and will 
present the fca /dependence of fco /for 6 = eB /= 
5.12 


1. Special case: T = 0, n = 0, 7 0 

According to the definitions of and 'Df~^ from 

(|IV.7I1 with from (IV.71) . the energy disper¬ 

sion relations of positively and negatively charged massive 
fermions are given by 

fc2(l - A+)2 - fc32(l - AXf - m2(i + c+)2 = 0, (V.8) 
and 

kl = kl + ml, (V.9) 


At zero temperature, the free parameter eB is scaled with the 
fermionic mass rriq instead of with the temperature T. 
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T=0, n=0, b=5 


a) T=0, q=+1, n=1, b=5 




k^/mq 


ks/mq 


FIG. 9: (color online). The ks/rug dependence of kojmg for 
positively charged (thick blue curve) and negatively charged 
(thin red curve) massive fermions in the LLL (n = 0) for b = 
eB/rrig = 5. They arise from the dispersion relations (IV.8II for 
q — +1 and (IV.911 for q = —1. 


respectively. 

In Fig. ini we have plotted the k-^jmq dependence of 
kq/mq for these fermions. As we have argued in Sec. lIVi 
two energy branches, appearing in this figure, correspond 
to positively charged fermions with positive spin (thick 
blue curve) and negatively charged fermions with nega¬ 
tive spin (thin red curve). Since fermions with negative 
charges receive no one-loop self-energy correction in the 
LLL, their dispersion relation is identical with the disper¬ 
sion relation of free fermions with kq/niq = 1 for k^ = 0 
(see the thin red curve in Fig. 0). On the other hand, as 
in the finite temperature case, the fact that the rest mass 
of positively charged fermions kq/niq > 1, is related with 
the dynamical creation of mass in the presence of strong 
magnetic fields, where the system is dominated by LLL. 
As it turns out from Fig. 0 for each fixed fcs/mq, the 
energy of positively charged fermions is larger than the 
energy of their negatively charged counterparts. This is 
also related with the same phenomenon of dynamical mass 
generation through the mechanism of magnetic catalysis 


2. Special case: T = 0, n = 1, m, / 0 

As aforementioned the one-loop corrected propagator 
of massive fermions in HLL is given in (IIV.8() - (IIV.11I) . 
To determine the energy dispersion relations of these 
fermions at zero temperature, we have to solve det(|f„ — 
niq — si^^) = 0 from (|IV.4I) . Alternatively, we can find 
the roots of the denominator Vn^ from pV.lOp with 
given in (IV.7() . Here, the denominator Vn '^, 
being a quartic function in , can be decomposed as 

(V.IO) 



ks/mq 


FIG. 10: (color online). The ks/mg dependence of ko/rrig for 
positively (panel a) and negatively (panel b) charged massive 
fermions for n = 1 and b — eB/m^ = 5. The results arise 
from the dispersion relations (IV.1511 corresponding to = 0 
(thick blue curves) and = 0 (thin red curves). 

with 

p(9) = (d(«)2 _ i = l,2, (V.ll) 

and 

4<?)2 ^ ^(q) _ l^^iq)2 ^ (y^^2) 

The functions and in (IV.121) are defined by 

0(9) = _ (a^9)2 _ 49)2 _ 4^)2^ |'^G)2 _ ^G)2 _ ^G)2^ ^ 

(V.13) 


with from (IV.7^ . The one-loop corrected 

fermion propagator in HLL is therefore given by 


SL^H~k) 


M^^\k) ^ M^^\k) 
v['‘\k) vi‘‘\k) ’ 


(V.14) 
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with Mn^ = ^Mn'^ and Nn'^ defined in 

(IIV.9I) . The energy dispersion relations, corresponding to 
the above two denominators, arise from = 0 , i = 1 , 2 , 
and read 

rf(9)2_^(9)2 ^ 

with from (jV.7ll and d!f\i = 1,2 from (IIV.12I1 . 

In Fig. nni we have numerically determined the fca jrfiq 
dependence of ko/mq, for n = 1 and b = eBlrriq = 5 as 
well as for g = +1 [Fig. [TOT al] and g = —1 [Fig. [TOT bl]. 
The thick blue and thin red curves in Fig. [TUI denote the 
solutions for = 0 and = 0, respectively. Here, 
in analogy to the hnite temperature case, described in Sec. 
Ill A[ two separate energy branches appear for positively 
and negatively charged particles. Their appearance is an 
indirect consequence of double spin degeneracy in HLL, 
in contrast to LLL, where the fermions have, depending 
on their charges, either positive or negative spins (see Fig. 

HI). 

Because of the similarity between the mechanism lead¬ 
ing to two energy branches in Fig. |U]at finite T and zero 
B, with the case discussed in the present section at zero 
T and nonzero B, the new fermionic excitations whose 
energy dispersions are demonstrated in Fig. [TUI are re¬ 
ferred to as cold magnetized plasminos. Let us notice at 
this stage, that, whereas hot magnetized plasminos ap¬ 
pear only in the limit of soft momenta ^3 <C T and weak 
magnetic fields eB <C (see Figs. [7|and[8|), as described 
in Sec. EKa), cold magnetized plasminos appear in the 
presence of moderate magnetic fields for all positive mo¬ 
menta fcs, and even in the massive case.^^ 


VI. SUMMARY AND CONCLUSIONS 

Plasminos are known to be collective excitations that 
appear in addition to normal fermionic modes in QED 
and QCD plasmas at finite temperature. Historically, 
they are shown to arise as one of the poles of the one- 
loop corrected fermion propagator at finite temperature 
T. The latter is shown to consist of two poles leading to 
two different energy dispersion relations. In the limit of 
vanishing fermionic mass and for small enough momenta 
of propagating fermions k, i.e. for |k| <C T, these energy 
dispersion relations belong to two fermionic modes with 
positive energies but opposite helicity to chirality ratios. 
The dynamical creation of additional collective modes in 
hot plasmas can be brought in close relation to the bro¬ 
ken Lorentz invariance, induced by the preferred reference 
frame defined by the heat bath. 


Moderate magnetic fields are necessary, because in a strong mag¬ 
netic field the fermionic system is solely dominated by the LLL, 
where, according to our previous arguments, no splitting occurs. 


A uniform magnetic field aligned in a fixed direction 
defines also a similar privileged reference frame. One of 
the consequences of such a frame is the appearance of 
certain anisotropies in the dynamics of fermions in the 
longitudinal and transverse directions with respect to the 
direction of the external magnetic field. Another conse¬ 
quence is the appearance of nontrivial collective modes in 
the spectrum of Dirac fermions in the presence of such 
a constant magnetic field. In the present paper, we ex¬ 
plored the quasi-particle spectrum of cold and hot QED 
plasma, by evaluating the general structure of the one- 
loop corrected propagator of magnetized fermions at zero 
and nonzero temperatures, and looking for its poles. 

To this purpose, we first determined in Sec. [HH the 
general structure of the free fermion propagator Sn'^ for 
positively (g = - 1 - 1 ) and negatively (g = — 1 ) charged 
fermions in the presence of a constant magnetic field B. 
In the momentum space, the free propagator is, in partic¬ 
ular, given in terms of the Ritus momentum fc„, where n 
labels the Landau levels [see (IH.41I) for free fermion prop¬ 
agator in the momentum space]. As it turns out, whereas 
HLL with n ^ 0 are characterized by a double spin degen¬ 
eracy, in the LLL with n = 0, the magnetized fermions 
have either positive or negative spins. This fact plays 
an important role in determining the correct spectrum of 
Dirac fermions in a constant magnetic field. 

In Sec. m we then determined the general structure 
of the one-loop self-energy of magnetized fermions, Tjn '^, 
at zero and nonzero temperatures. The main results are 
presented in (IHL18I) for T = 0 and in (IHL32|) for T ^0. 
A comparison between these two expressions shows that 
the isotropy between the components of the longitudi¬ 
nal momentum ky = (fco,0,0, fca) at zero temperature is 
removed at finite temperature. This is because of the 
aforementioned breaking of Lorentz invariance at finite 
T. We computed Tjn'^ for g = ±1 in terms of certain co¬ 
efficients, and , and determined these 

coefficients separately for T = 0 and T 7 ^ 0 up to a num¬ 
ber of integrations and a summation over Landau levels. 
To explore the Dirac spectrum for soft momenta ^3 <C T 
and in the limit of weak magnetic fields eB <C , we also 
evaluated the above coefficients in a HTL expansion in a 
weak magnetic field. 

The one-loop corrected propagator of magnetized 
fermions Sn'^ is then determined by combining the free 
fermion propagator Sn'^ and the one-loop fermion self¬ 
energy in the standard way [see Sec. [IV]. Here, 
we have not distinguished between the zero and the fi¬ 
nite temperature case. The main results are presented in 
(IIV.6I) - (IIV.7I) for n = 0 and in (IIV.8|) - (IIV.11I) for n 7 ^ 0. 
We considered two special cases (n = 0,mq 7 ^ 0) as well 
as (n ^ 0,rnq = 0 ), and determined the spectrum of 
fermionic modes by computing the eigenvectors of the nu¬ 
merators of Sn^. We showed that they are, in particular, 
eigenvectors of E 3 = diag(cr 3 , 0 - 3 ), defined by the third 
Pauli matrix (T 3 . Hence, in the LLL, the spectrum con¬ 
sists of positively (negatively) charged fermions with pos¬ 
itive (negative) spins, while in HLL the expected double 
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spin degeneracy occurs. In addition, massless fermions 
with n ^ 0 turned out to have well-defined left (negative) 
or right (positive) chiralities, as is described in (IIV.26I) . 
but they are not eigenstates of the helicity operator, in 
contrast to the standard T 0 and B = 0 case. In 
the massless case, the general structure of in HLL is 
presented in (IIV.20|) - (IIV.22l) .The appearance of two de¬ 
nominators X>L and Vji in (IIV.20I) leads to two different 
energy branches, similar to the case of T ^ 0 and B — 0 
[see (III.121) ]. although the corresponding collective modes 
have different properties. 

In Sec. El we numerically determined the energy dis¬ 
persion relations of these magnetized collective modes at 
zero and nonzero temperatures. In Sec. IV Al we con¬ 
sidered first the case of T ^ 0 and B 0, and de¬ 
termined the fca/T dependence of ko/T in two different 
regimes of k^/T < 0.2 (soft momenta) and k^/T > 0.2 
(hard momenta). To solve the energy dispersion relations 
arising from the poles of the corresponding propagators, 
we numerically determined the above mentioned coeffi¬ 
cients and in these two regimes, for 

fixed values of eBjT^ and for massive as well as massless 
fermions. In the LLL, only one energy branch arises in the 
whole regime of k^/T. It belongs to positively or nega¬ 
tively charged fermions with positive or negative spins. In 
contrast, according to the results from Fig. [71 additional 
excitations, referred to as hot magnetized plasminos, ap¬ 
pear in HLL, where left- and right-handed fermions have 
both positive and negative spins. In the limit of weak 
magnetic fields eB <C and soft momenta k^ <C T, 
these collective modes have different energies, while for 
larger values of k^/T they seem to have the same energy 
dispersions. This is in contrast to the results arising for 
nonzero B and zero T, discussed in Sec. EH In this case, 
the results for n = 0 and rriq ^ 0, are qualitatively the 
same as in the finite temperature case. For n ^ 0, how¬ 
ever, cold magnetized plasminos appear in the presence 
of moderate magnetic fields and for all positive momenta 
ks, even in the massive case [see Fig. dU]. 

Let us finally notice, that the group velocities of prop¬ 
agating collective modes can be determined from their 
energy dispersion relations for (T ^ 0,B ^ 0), and might 
have applications, e.g., in the physics of heavy-ion colli¬ 
sions. At finite T and zero B, the appearance of a mini¬ 
mum at some finite value of the momentum in the energy 
dispersion of plasminos, leads to a vanishing group veloc¬ 
ity for the collective modes. The latter has been inter¬ 
preted as the appearance of Van Hove singularities 23|, 
e.g. in the low mass dilepton production rate in the QCD 
plasma . The sharp structures arising in this quantity 
are known to provide a unique signature for the presence 
of deconfined collective quarks in the quark-gluon plasma 
M- It would be interesting to determine the dilepton 
production rate at high temperature and in the presence 
of moderate background magnetic fields. These are be¬ 
lieved to be produced in early stages of heavy-ion col¬ 
lisions [13 , and, because of certain medium effects are 
assumed to be approximately time-independent (^ . We 
will postpone these kind of phenomenological studies to 


our future publications. 
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Appendix A: A useful formula 

In this appendix, we will compute the integral 

^+oo 


in.P — 


/ H-OO ^ 

dxie~^p'^^ f+{xi,k2)fi{xi,p2)- (A.l) 

-OO 


Here, according to (III.291) . 

/+(a;i, *2) = a„e Hn 


'Xi - ilk2 


(A.2) 

with iq = \qeB\~^/‘^ and an = (2”n!y^£g)“^/^. To this 
purpose, we use the following representation of the Her- 
mite polynomial Hn(x) 


n! f dt 




tn+l 




arising directly from 


-A+2tz ^ 




(A.3) 


(A.4) 


Plugging (IA.3I) in O, we arrive first at 
a„a£nW. / dt du 

In p = 


(27ri)^ J 

^ g-(P+u'’) g-2^, (tfe2+up2) g- 

X / dxie 1 ^ . (A.5) 

J —00 

The integration over xi can be performed by quadratic 
completing the square, and performing the resulting 
Gaussian integration over xi. This results in 

_ -f(fc2-P2)2-ep'?]^ -ipi(/c2-eP2)^ 

“ (2^i)2 ® 

X J„£(fc2,P2,p'l), (A.6) 


with 


J f = J) ^'2tu+tla(p2-k2-ip',)+uia(k2-p2.-ip'-,) 

(A.7) 

To determine Jni, we will first perform the integral first 
over t and then over u. Then, integrating first over u and 
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then over t, and comparing the two resulting expressions, 
we will eventually arrive at the generalized formula for 
Jni{k2,P2,p'i)- 

The integration over t in (IA.7I) can be performed, using 
the Cauchy formula 


dt . 27ri dP . 


2=0 


which leads to 


/ 


tA(u) 


n\ 


(A.8) 


(A.9) 


Setting A(it) = 2M+£q(p2—fc2—*K), plugging the resulting 
expression in (IA.7I1 . using the binomial series identity 


( 2 m + iq{p2 - k2- ip'i)y 


= E 


r=0 


{2uniq{p2 - fc2 - tp[)r-^ 

and eventually replacing -^ith 

^u{k,-p,-^p[) _ ^ U^ik2 -P2- ip[y 


(A.IO) 


i=o 


J! 


(A.ll) 


we arrive first at 


. n 00 ryr nn—r-\-i 
2Tn 2 tq 


J 


r—0 j—0 

X {P2 -k2- -P2- ip'iY 

du 


^e-r-j+l■ 

Integrating then over u, 
du 


u‘ 


l-r-j+l 


— 27rzJr+j,t, 


(A.12) 


(A.13) 


we obtain 


-k2- ip'iT{k2 -P2- ip'iY 


with 


n\£\ 

xKne{k2,p2,p'i), 


Vi (z), 


(A.14) 

(A.15) 


and z = -^[{k 2 — P 2 Y +^1]- Here, W^(z) is the confluent 
hypergeometric function of the second kind [l^, defined 

by 


Let us now consider (IA.7I) , and check what would happen 
if we computed Jni by integrating first over u and then 
over t. In the above case, where the integration over u is 
performed after the integration over t, we have, according 
to (IA.I3|) . £ = r+j. For j > 0, arising from (IA.12I) . we get 
therefore £ > r. On the other hand, it is clear from (IA.12I) 
that r < n. Combining these results, we arrive at n < £. 
This fixes the upper limit in the summation over r, which 
is given by n equal to min(n,^). If we computed Jni by 
integrating first over u and then over t, we would arrive at 
a summation over r from r = 0 to r = £, and, according 
to the above argument, £ would be equal to min(n,^). In 
other words, the most general expression for is given 
by (IA.14I) with n and £ replaced by to = min(n, t) and 
M = max(n,^). Plugging now this expression for Jni in 
(IA.6I) . the most general expression for Int reads 

In,l = 

X {P2 -k2- ip[)^{k2 -P2- ip'x)^ 
xe-fz—(A.17) 

where Ani = (2"+^n!£!)“^/^, z = ^[(^2 - P 2 Y +p'?], 
TO = min(n,£) and M = max(n,I'). In Sec. HID the 
above relation will be used to evaluate the integration 
over xi and yi in (IIIL6I) with k 2 ,P 2 ) from 

(|III.10I) . Here, we present only the result for one typical 
combination, 

r rt _ _ ^ _ —z M—m\u—m (~\\2, 

7j!(£ _ 1). ^ l^M—m+lw/J ’ 

(A.18) 


C 9 2 

with z = -^{{k 2 — P 2 ) +^1]; ur = min{n,£ — 1) and 
M = max(n,£ — 1). All the other integrals appearing in 
(|HI.6I) with from (IHL7|) are performed in the same 
way. 
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